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Abstract 

The kinetics of the diffusion-controlled chemical reactions A+A+...+A — > that occur 
at catalytic centers periodically arranged along a straight line is considered. Modes of the 
behavior of reaction probability W(t) were studied at different times and different densities 
of the catalyst. Within the Smoluchowski approximation, it was rigorously proved that 
at large times the function W(t) is independent of the lattice period. This means that, 
in the given asymptotic mode, the probability of the reaction on a lattice with a catalyst 
placed in each lattice site is the same as on a lattice with a catalyst placed in sparse sites 

In this work, we discuss the kinetics of the polymolecular irreversible diffusion-controlled 
chemical reactions A + A + ... + A^0 that occur at catalytic centers arranged in a periodic 
lattice in space. To our knowledge, the problem of the decrease in the amount of the reactant 
with time in catalytic diffusion-controlled chemical reactions A + A—>-0 was first studied within 
q ' the random phase approximation in [|TJ . 

The problem considered in this article "rests on two whales": on the one hand, it has 
.£h ! characteristic features of the dynamics of diffusion-controlled chemical reactions, and, on the 
other hand, it is typical for the adsorption kinetics of particles diffusing in the medium of 
stationary "traps" (absorbers). Before going into the description of the model under study, 
let us briefly formulate the main features of its fundamental problems: diffusion-controlled 
chemical reactions and kinetics of the adsorption of particles on traps. 

For several decades, bimolecular diffusion-controlled reactions A + A — > B served as the 
basic model for studying the kinetics of chemical reactions on the microscopic scale. The main 
ideas of the description of chemical reactions in the medium of diffusing reactants were formu- 
lated by Smoluchowski in 1917 0. Mean field methods, which were proposed by Smoluchowski 
and further developed by his numerous disciples (see, e.g., review are widely used up to now 
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and yield reasonable results for problems related to the determination of the rate of chemical 
reactions and change in the concentration of the reactant with time. However, in the case of 
strong stirring, i.e., when the diffusion of particles is significant, mean field methods become 
inappropriate because, along with pair interactions of particles, correlations that occur at triple, 
quadruple, etc. particle collisions should be taken into account. In works that date back to the 
1980s, fluctuation effects in the kinetics of chemical reactions were included Later the range 
of problems under study was signifi- cantly expanded by the inclusion of problems of chem- ical 
kinetics involving macromolecules, so-called transport-controlled chemical reactions P|. How- 
ever, the comparison of the results that were obtained within the Smoluchowski approximation 
with exact solutions for some model systems, e.g., for a one-dimensional bimolecular reaction 
A + A — > (without catalyst) demonstrates that the mean field method provides cor- rect 
scaling for the probability of particle decay at large times || which differs from the exact 
solution only by a numerical coefficient. 

The problems of determining the evolution of survival or decay probability for particles 
diffusing in d-dimensional space in the presence of stationary randomly arranged absorber traps 
is the subject of many studies in both the physicochemical and mathematical literature that 
have appeared over the last two decades or more. The great interest in this topic is obviously 
caused by its significance in studying such physical chemistry problems as photoconductivity 
and photo- synthesis, the kinetics of binding biopolymers with ligands, and the adsorption 
of polymer molecules on surfaces with chemically active regions and on colloidal particles in 
solutions and gels. Works dating back to the 1970s 0] demonstrated that the probability of 
the survival of independent diffusing particles in a medium with absorbing traps at large times 
cannot be described within the Smoluchowski mean field approximation and is controlled by 
the diffusion time of a particle in a cavity of typical size free from traps. The method that 
provides correct asymptotic relationships is identical to the optimal fluctuation method by 
I.M. Lifshitz for estimating the density of states of a disordered system in the vicinity of the 
band edge (so-called Lifshitz singularities). Thus, works |7|, [3], f| stimulated the interest of 
mathematicians in the problem of particle decay in a medium of randomly arranged absorbers. 
Subsequent studies of model systems demonstrated the relation of this problem to percolation, 
the statistics of polymers in a medium of random obstacles, and some aspects of supersymmetric 
quantum mechanics P, 3; PI; CT- 



Now, with the general notion as to what type of problem of statistical physics our problem 
belongs, we can turn to the discussion of the model, which, as noted above, was formulated in 
[p]]. Let us consider a region ("reservoir") in ^-dimensional space containing an ensemble of 
identical particles A of finite size involved in Brownian motion and some quantity of stationary 
catalytic centers (traps). In the course of a random walk, there is the probability of encountering 
of m particles. If this encounter takes place at a trap, the elementary event of the chemical 
reaction 



A + A + ... + A^0 

occurs and all the particles decay, i.e., are removed from the reservoir; if particles encounter 
outside the trap, the reaction does not occur. The task consists in calculating the probabil- 
ity W(t\m) of the particle decay within the given time interval t depending on the spatial 
arrangement of traps. 
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The probability of particle decay W(t\m = 2) for a bimolecular chemical reaction was 
estimated in |J within the Smoluchowski mean field theory for the periodic and random ar- 
rangements of traps, and the effect of catalytic centers was included in || within the random 
phase approximation. In our work, we restrict our consideration to the case of d — 1 and 
the periodic lattice of traps. For this model, the distribution function W(t\m) can be deter- 
mined exactly without any additional assumptions. We emphasize that by the exact solution 
we mean only the possibility of the rigorous inclusion of the effect of the spatially distributed 
catalyst on the probability of particle decay W(t\m = 2). Note that, from the standpoint of 
real physical problem when a finite concentration c of a reacting particles occurs in the system, 
we still remain within the Smoluchowski approximation. At the end of the article, we discuss 
the range of applicability of the given mean field approach. Although one- dimensional, the 
problem in hand has a particular physicochemical application. Chemical reactions that occur 
at cathodes in the presence of a catalyst (platinum particles) can be readily described within 



the model under consideration [II], [I2| . This is related to the fact that on the microscopic 
scale the cathode is a porous structure with one-dimensional channels in which the catalyst 
is deposited and the diffusion-controlled chemical reaction occurs. A problem of fundamental 
importance consists in the optimization of the adsorption process: increasing the rate of the 
chemical reaction and decreasing the concentration of the expensive catalyst. Let L be the 
distance between the adjacent traps in a straight line. We denote the time dependence of the 
coordinates of the particles under consideration by X\(t),X2{t), ,x m (t). The condition that all 
particles at the instant of time t occur at an arbitrary trap with the coordinate nL is as follows: 
Xi(t) = x 2 {t) = = x m (t) = nL . It is easy to realize that the random walk of m independent 
particles is equivalent to the effective one-particle random walk in the m -dimensional Euclidean 
space £ m = (xx, x 2 , , x m ). Thus, the probability W{t\m) of the decay of m identical independent 
particles randomly wandering along a straight line at their simultaneous encounter at any of 
the traps is equivalent to the probability of the first occurrence of the random walk in the S m 
space at the instance of time t in an arbitrary trap located in the straight line defined by the 
equation x\ = x 2 = = x m . In our case, this statement of the problem is a starting point. The 
recurrence equation for the probability Wt(x) = W(t\m, x) to find the random walk over lattice 
edges in the m-dimensional space S m in the point x at the instance of time t has the form 

iy m (x) = D W t (x + u)77(x + u), (1) 

u 

where D = is the diffusion coefficient in m-dimensional space, summation is performed over 
the nearest neighbors, and n() is the discrete ^-function, which is defined as follows: 

, x _ J 0, if x coincides with the trap; , . 

| 1, in all the other cases. 

As usual in translation-invariant problems, it is convenient to transform to the momentum 
Fourier representation: 

^( k ) = ^E^) eikx > ( 3 ) 

X 

where 

Wi(x)= Yl ^(k)e" Jkx , *i = ±^T (4) 

|fel,...,fc m |<7T 
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As a result of the Fourier transform, Eq. ([!]) yields 

= h E e " iku E ^(q)^( k - q), (5) 



2m 

U q 

where 

ry(k) = 5(k) - ^ ( fi ) 

|nL|<7V/2 

In Eq.(^|), we assumed that the straight line in which the traps are located coincides with the 
[0, %i) axis in the E m space. Let us use the equation 

\ 

2L 



N i i . L \ / 



in which, for the sake of defmiteness, we assume 

IM -= 2L- Usill S tlle l£Ltter 

expression, 

Eq.@ can be written in the form 



W t + 1 {k) = -{cosk l + ... + cosk m )\w t {k)- — I— V W t (k 



I - -t-,Q2, -Am ) }■ (7) 



\n\<L/2 

where 

-7T < fc m } < 7T, -7T < {g 2 , g m } < 7T 

Equations of this type are conveniently solved by the generating function method: 



W '(k, s )=f>(k) s ', H'W.l/^. (8) 



t=o c 



Multiplying both parts of expression (0) by s*, after simple algebraic transformations we obtain 

£/ H/p(k) 



where the following designations are introduced: 

A(k,s) = — (cos k\ + ... + cos k m ), (10) 

777 



12 8m \ 

\n\<L/2 



2im 

92) • • • j Qn 



Let us select uniform starting conditions, i.e., Wo(k) = 6(ki) . . . S(k m ) and proceed to limit 
N -> oo: 



gk) A(M i / v / /• d Q2 . . . 



-1 



1-A(k,s) 1 - A(k, s) (1 - s) \^ J J l-^(cosgi + ... + cosg„ 



(11) 
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where q\ takes on the values q\ = k\ j 1 , |n| < L/2. The second term, which enters into 

the expression with the minus sign, describes particle decay at traps. We are interested in the 
zeroth harmonic of this term, which determines the probability of the reaction 

w(s) = irh? x (E / ■ ■ ■ / 1 - ^(cret + d " m + c°s 9m ) ) ' • < 12 > 

For further consideration of the problem, we must evaluate the integral 

ti \ f f dq 2 ...dq m 

I W= "' 1 7 T V' 13 

J J a - (cos q 2 + ■ ■ ■ + cos q m ) 

where a = — — cosgi. This function can have poles only on two points a = ±(m — 1), 
i.e. si } 2 = cos 9l -^ra,!) ; $1,2 > 1- Integral (|13|) in these regions is determined by the values 
of qi, which are nearly zero; therefore, cosQ in the integrand can be expanded. The change 
cos qi ~ 1 — corresponds to the changing from the lattice to the continuous limit with 
respect to the ith coordinate. In the vicinity of the point a = m — 1, e — a — m + lwe have 



Changing to spherical coordinates (S m i is the area of the sphere in the (ml)-dimensional space, 
A is some constant, A ~ ), we obtain 



A q m ~ 2 dq 
2e + q 2 



2^-1/ hr^- (15) 



The values of I m (e) at different values of m are 

27T 2 



— m = 2- 

87rlog(l + f), m = 3; (16) 
~ A™- 3 , m > 4. 

The probability of decay at traps within these designations is 

777 

W (s) = ^ • (IT) 

p=-L/2 

The integrand in the function W(t) has singularities s* and (1 — s) 2 . We are interested in the 
behavior of W(t) governed by the vicinity of the point s = 1, which corresponds toa = m— cos q\. 
Hence, the expression for the generating function W(s) within a numerical multiplier can be 
rewritten as follows: 

W(s) — . (18) 

(1- S ) 2 E Im(f -cos^-m + l) 

p=-L/2 
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The time dependence W(t) is restored from the generating function W(s) by inverse Laplace 
transform (§). The contour integral depends on singular points of the integrand, and the pole 
at the point s = 1 makes the largest contribution. Recall that the probability W(t) of particle 
decay unambiguously characterizes the decrease in the reactant concentration with time and, 
thus, determines the effective rate constant of the chemical reaction. Further calculations will 
be based on the Tauberian theorem, which readily provides asymptotic estimations in the cases 
of interest without explicit use of the inverse Laplace transform. 



Reactions with m = 1 and m > 4. Retaining only the divergent part in expression (|12|) , we 
obtain the asymptotics 

W m=1 (t -> oo) ~ y/i. (19) 



In the case of m > 4 the result is also readily obtained because I m ~ const in the region 
s — 1. Therefore, t 
which corresponds to 



of s — 1. Therefore, the behavior of the function at large times is governed by the pole n~p ? 



W m > 4 (t -> oo) ~ t. (20) 

Reactions with m = 2. The behavior of W(t) at m = 2 is more interesting. In this case, the 
generating function W(s) of the absorption probability is as follows: 

W(s) — . (21) 

(l_ s)2 E ( 2 -scos 2 -f-s)- 1/2 

p=-L/2 

Applying the Tauberian theorem and replacing the sum with the integral, we obtain 

Wm - 2it) ~ VU^Tf [log4-log( i + ff )]' (22) 



where a — y + \- In the two limiting cases 



W m=2 (t) ~ ( ~ ^ f L2 (23) 

[ LlogSt' 1 ^ 1 ^ 7T 2 - 

Reactions with m = 3. Expression ([T^) for I m= s gives the correct relationship; however, it 
contains an undefined constant, which complicates estimations. By the integration of Eq. fljjD 
without expanding cosines in series, we can obtain the more accurate formula 

1(a) » vr log ( - - 16 ^ . (24) 



3 — 2s — s cosg / 

The asymptotics for t — > oo can be obtained directly by applying the Tauberian theorem to 



generating function (17), which leads to the following expression: 



C + log*-log(f + 
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where C = 5L log 2 + 4 log ) - 4 - log f . 

Absorption probabilities that were obtained above are insufficient for the exact calculation 
of the reaction rate but are directly related to the Smoluchowski constant Ks mo i(t) = 
Indeed, this constant by defi- nition is the probability of the decay of a solitary particle (i.e., 
without inclusion of cooperative effects) at the catalyst trap. Within the Smoluchowski approxi- 
mation, the concentration of particles in the polymolecular reaction of m particles is determined 
from the kinetic equation 

^ = -K Sm oi(t)C tr C m (t) (26) 

After simple calculations, we obtain (Co = C(t = 0), m > 1) 

C(t) = (Co 1 -" 1 + (m - l)C tr W{t)) - 1/(m ~ 1} . (27) 

The Smoluchowski approximation is valid in cases when stirring in the system due to the 
diffusion of particles is a slower process than the chemical reaction event. That is, if the motion 
of particles is nearly absent, the chemical reaction event involves only the particles that are 
occasionally the closest to each other. This is the reason why the Smoluchowski approximation 
is adequate at a low density of particles in the system, when the free diffusion path of particles 
is sufficiently large compared to the time of the chemical reaction event. To demonstrate the 
validity of the mean field approximation in the problem under consideration for any distance L 
between the catalytic centers, it is sufficient to show that the Smoluchowski theory works in the 
two limiting cases: (a) at L — 1, i.e., when the reaction occurs each time on the encounter of a 
pair of particles and, consequently, is independent on the position of the catalytic center; and 
(b) at L —>■ oo, i.e., when the system contains a solitary catalytic center. The exact solution of 
the many-particle problem in case (a) for bimolecular reactions was reported in [[| ^3J. As noted 
in the introduction, the exact and mean-field solutions have the same asymptotics and differ 
only by a numerical coefficient. For the case (b), the decrease in concentration c(t) because of 
the chemical reaction can be easily estimated (m > 2): 

dCchemjt) . . . . C ,-l/( TO -l) / oq n 

— = -K chem c chem (t) } =^> c chem (t) ~ — t (28) 

dt J\ chem 

In order to find the decrease in concentration c(t)dis because of diffusion in the case of the 
reaction of m particles at a center, we must solve the diffusion equation in m-dimensional space, 
where each coordinate corresponds to the concentration of one of the particles involved in the 
reaction. As can be easily seen, 

Cdiff(t) ~ c (m = 2) ^ 



Cchem(t) ~ c const (to = 3, 4, 



For all m > 2, the following condition is fulfilled: after some instant of time t = t(co, K c h em ), 
the concentration because of diffusion is larger than the concentration because of the chemical 
reaction, and, consequently, the rate of diffusion is lower than the rate of the chemical reaction, 
which means that the Smoluchowski approximation is valid. 
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Of particular interest are the relationships W m (t) at large times ( t — > oo), namely, expres- 
sions (0), (p3|) and (p5|). As can be seen from the corresponding formulas, the probability 
of reaction W m (t) is independent of the period of the lattice of catalytic centers. This means 
that in the given asymptotic mode the probability of the reaction on the lattice with the cat- 
alyst placed in each lattice site is the same as on the lattice with the catalyst placed in sparse 
sites. Recall that this result was first formulated in |IJ within the random phase approximation. 
Thus, our work can be considered as rigorous proof of the effect, which has promising technical 
applications. We are grateful to G. Oshanin for helpful discussions of the work. The work was 
supported by the Russian Foundation for Basic Research, project No. 00-15-99302. 
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